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Abstract 



We attempt the numerical construction of an effective action in three dimen- 
sions for Ising spins which represent the Wilson hues in the four-dimensional 
SU(3) gauge theory at finite temperature. For each configuration of the gauge 
theory, each spin is determined by averaging the Wilson lines over a small 
neighborhood and then projecting the average to ±1 according to whether 
the neighborhood is ordered or disordered. The effective Ising action, deter- 
mined via the lattice Schwinger-Dyson equations, contains even (two-spin) 
and odd (one- and three-spin) terms with short range. We find that the 
truncation to Ising degrees of freedom produces an effective action which is 
discontinuous across the gauge theory's phase transition. This discontinuity 
may disappear if the effective action is made more elaborate. 
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I. INTRODUCTION 



In studying the fluctuations of a complex statistical field theory, it is frequently useful to 
define simple, effective degrees of freedom. A wise choice of these degrees of freedom allows 
one to focus on specific physics. We are interested in the first-order phase transition of 
the SU(3) gauge theory specifically in surface phenomena such as the surface tension 
between the phases and the stability of bubbles [^]. Since the transition is first-order, 
all correlation lengths stay finite on either side. This means that there should be no delicate 
issues related to renormalization-group fixed points and the symmetries that characterize 
them. We are thus led to degrees of freedom which simply specify whether a given region is in 
the confined or unconfined phase. Assigning values of ±1 according to the two possibilities, 
we are led to an effective theory of Ising spins. 

This exercise is familiar from study of the liquid-gas transition [0. One assigns values 
to a local Ising spin a according to the local density p of the fluid. At the first-order phase 
boundary T = T* the density is discontinuous, and this is reflected in a discontinuity in 
the magnetization (cr) of the effective Ising theory. The simplest Ising action has only one 
term that is odd in the spins, namely, the magnetic field term ha. The first-order transition 
occurs perforce at h = 0, and one identifies the liquid-gas phase boundary with the segment 
of the T axis between the origin and the Ising critical point T^^^^^. As one slides along the 
phase boundary towards the liquid-gas critical point, the discontinuity in p decreases until 
it reaches zero at the critical point; correspondingly, the discontinuity in (a) vanishes as one 
approaches Tjf^^. 

The simplest Ising model, however, is incapable of describing an arbitrary first-order 
phase transition. An easy way to see this is to note that its Z{2) symmetry implies that 
{o-)h^o+ = —{o')h^Q-. Since the correspondence between a and p depends on an arbitrary 
assignment in the first place, there is no reason to suppose that this symmetry reflects 
accurately the values of p on either side of T*. Hence there must be more than one odd term 
in the Ising action, in order to move the transition away from h = and thus to destroy the 
Z{2) symmetry about the transition. 

For the SU(3) confinement transition, we seek to characterize a neighborhood of a site n 
as confining or non- confining. We naturally settle on the Wilson line as the quantity 
which does this. As an order parameter of the Z[3) symmetry of the Euclidean theory, L 
is discontinuous at the transition, with (L) = in the confining phase below the transition 
and (L) 7^ in the plasma phase above. We use the modulus of L^, suitably smeared, to 
assign a value to = ±1. (This smearing reflects the fact that confinement is a property 
of a neighborhood, not of a point.) Running a Monte Carlo simulation of the gauge theory, 
we generate configurations of which translate into configurations of a^. We then use 
the Ising model's Schwinger-Dyson equations to determine an approximation to the 
effective Ising action S'efT[cr]. 

This definition of a as a function of L integrates over the Z{3) dynamics of the Euclidean 
theory. Since we are interested in understanding bubbles of the confining phase in the plasma 
and vice-versa, distinguishing among the three ordered phases is unnecessary. Previous work 
IP has taken the opposite approach, projecting the complex Wilson line onto 3-state Potts 
spins r = exp27rnz/3. This makes it harder to identify bubbles of the disordered phase 
because they will only be visible when calculating averages of r over sizable neighborhoods. 
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Since our cr spins depend on the magnitude of L, they show directly those places where \L\ 
is small and hence disordered. 

Since Z{3) domain structure plays a role in the confinement phase transition, one might 
be concerned that domain walls should somehow be represented in the effective action. We 
note that calculations ^ have shown that, near the transition, the disordered phase wets the 
ordered phases, and moreover that two ordered domains with different Z{3) orientations will 
sandwich a disordered domain between them. We note also that the physics of fluctuations 
among the three ordered phases will influence the Ising couplings because the effective 
action comes from an integration over all other degrees of freedom in the gauge theory. 
Nevertheless, integrating thus over the Z{3) fluctuations may impair the ability of the Ising 
theory to represent the phase transition correctly. Our effective Ising action turns out to 
be discontinuous at the phase transition. Although this might be due to the small number 
of terms we permit in the action, it might stem from the use of Ising variables itself. We 
discuss this further below. 



II. DEFINING THE EFFECTIVE ISING THEORY 

We simulate the SU(3) gauge theory at flnite temperature by using, as usual, a Euclidean 
lattice of A^^ x A^^ sites, with the physical temperature given by T = {Nta)~^ in terms of 
the lattice spacing a. All results presented in this paper were obtained from lattices with 
Aj = 2. The gauge theory is governed by the Wilson plaquette action, 

= /? E Tr KU:^f,Ui:i,U:^ , (1) 

rt 

with Uf^ G SU(3). The order parameter for the conflnement phase transition is the Wilson 
line, deflned on a site n of a three-dimensional lattice via 

Nt 

Ln = Tr n ^(n,no) " (2) 
no=l 

As discussed in the Introduction, we will use L to deflne the effective Ising spins a which 
will label a neighborhood on the lattice as conflning or non-conflning. A flrst attempt might 
be to deflne 

The problem with this is that Ln fluctuates violently from site to site. Even deep in the 
connflning phase, where (L) = 0, is by no means conflned to a region around zero, and 
in fact fllls the entire wedge available to it in the complex plane (see Fig. |I]). a as defined 
by thus does not offer a good definition of a domain in the confining phase. 

The fiuctuations in Ln are reduced if it is averaged over a small volume. We define 
to be the average of L over the m x m x m block surrounding^ n. A glance at Fig. |l| 



^If m is even then n is a site of the dual lattice. 
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shows that L[^l discriminates well, on a local basis, between domains that resemble the two 
respective phases. Ll^^l , on the other hand, fluctuates too little about the volume average of 
L, so that, with a reasonable value chosen for To-, the a spins would lose all information about 
fluctuations and remain entirely ordered with a = ±1. We thus choose for insertion in 
(0) to calculate the a conflgurations; as shown in Fig. ^, we set = 0.8. 

With the deflnition of a in hand, we turn to the determination of the effective Ising action. 
In principle the action has an inflnite number of terms; we truncate it to a combination of 
a magnetic fleld term and two- and three- spin terms with range 2, 

S.s[a]=J2(^c.O\ (4) 

a 

where the seven operators 0" are listed in Table |I[ The two-spin operators and O^, as 
well as the three-spin operator O^, couple spins within distance ^/2■, the remaining operators 
reach out to distance 2. 

A Schwinger-Dyson equation of the Ising theory is derived by flipping a spin in the 
sum deflning the expectation value of some operator. For the operators in Table |, we have 

{dz) = -{dZexp2S^) (5) 
where we have deflned (9^ to be those terms in 0° that contain cxn, and 

5„ = E/5«(^S (6) 

a 

is the part of the action that contains a^- These are seven equations for determining the 
seven unknowns f3a- After generating an ensemble of a conflgurations via Monte Carlo 
simulation of the gauge theory, we determine Pa iteratively as solutions of (|]). 
As a consistency check, one may use the vacuum equation 

l = (exp25n) (7) 

or the Schwinger-Dyson equation for any other operator in the theory. A more satisfying 
check, however, is to run a direct Monte Carlo simulation of the Ising model with action 
to see if the expectation values of as computed in the gauge theory are reproduced. This 
was the procedure we followed. We calculated error estimates by subdividing the ensemble. 



III. RESULTS AND DISCUSSION 

We simulated the SU(3) gauge theory on a lattice of volume 2 x 16^. The conflnement 
phase transition is in the neighborhood [|TU| of /3 = 5.09, and we settled on the value 
(3 = 5.091 after seeing no tunneling between the coexisting phases in moderately long runs 
at that coupling. 

Straightforward application of the method described above gives an Ising action for a 
three-dimensional lattice of volume 16^. We show the couplings for this action, derived from 
ordered and disordered runs at /? = 5.091, in Table |T[ In both cases, the action contains 
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couplings with range a/3 and 2 (O^ and 0^, respectively) which are as strong as the shorter- 
ranged two-spin couplings [O"^ and O^) and compete with them in sign.0 This raises the 
suspicion that a longer-ranged Ising action is needed to reproduce the Ising configurations 
correctly, and that the range-2 action is too crude a truncation. This suspicion is confirmed 
by simulating directly the Ising model with the couplings just derived. As seen in Table fT|, 
comparison of (O") with the averages from the gauge configurations shows poor agreement.^ 

This problem was encountered by Deckert et al. |]lT| in a calculation of the effective action 
for the Z{2) gauge theory. A solution, noted in P], is to perform a block-spin transformation 
on the spins, so that the effective action has twice the range. We do this simply by decimating 
the Ising spins, already defined via smeared averages, to an 8^ sublattice. Solving the 
Schwinger-Dyson equations with the decimated configurations gives the couplings shown in 
Table The longer-ranged two-spin couplings, [3q and /?7, are negligible, as is the straight 
three-spin coupling (3^. Moreover, comparison of an Ising Monte Carlo simulation with the 
gauge theory (see Table |V|) now gives satisfactory agreement. 

The effective Ising couplings shown in Table |^ vary smoothly as (3 is varied on either 
side of the transition, but they are discontinuous across the transition. At /? = 5.091 we 
have, then, two actions, S'coid and S'hot, which are the limits of S'efr[cr] from the ordered 
and disordered sides of the transition. Curiously, we find that S'hot is at a point of phase 
coexistence, that is, at a phase transition between phases with (a) < and {a) > 0. We 
show in Table the expectation values of the seven operators (9° for S'coid and for both 
phases of Shot- For Shot, the phase with (cr) < describes well the expectation values in the 
gauge theory on the disordered side of its transition. The other phase of S'hot, of course, does 
not; neither does it describe the ordered phase of the gauge theory. The action S'hot thus 
"knows" that it describes a phase transition, but it is capable of describing correctly only 
one of the phases. S'coid describes the ordered phase well, and shows no phase coexistence. 

The discontinuity in the effective action is an example of the singularities that can result 
from renormalization group transformations. Griffiths and Pearce |jl2| noted that a blocked 
action might be a singular function of the unblocked couplings even though the blocking 
transformation is local. Later work found discontinuities in the blocked action associated 
with first-order phase transitions in the original action. It was conjectured that there may be 
different renormalizat ion-group flows resulting from the various metastable phases at a flxed 
coupling. In view of theorems proven by van Enter, Fernandez, and Sokal [|I^], however, 
such discontinuities are impossible in an effective action which possesses flnite range in the 
inflnite- volume limit. Our effective action, however, is approximate in that it contains a 



^Note that the magnetic field is h = —Pi, and that negative values for P2, Ps, Pe, and P7 indicate 
ferromagnetic couplings. 

^The violent disagreement for the ordered phase, including even the sign of the magnetization, 
suggests that the gauge theory's operator averages are to be sought in a metastable phase of the 
Ising action. We did not succeed, however, in reaching this phase with our Monte Carlo. In the 
disordered phase, the positive magnetic field h = — /3i prefers a positive magnetization, but the 
positive three-spin couplings P4 and (3^ compete with it and turn the magnetization negative. 
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small number of couplings. Adding longer-ranged and multi-spin terms to the action will 
bring consistency with the theorem of van Enter et al. in one of two ways: Either the 
couplings will become continuous |T^, or the effective action(s) will acquire too many non- 
local terms in the infinite- volume limit, meaning the statistical measure is non-Gibbsian. In 
the first case, we will have an effective action well suited to describing the phase transition; 
in the second case, the conclusion will be that an Ising description of the phase transition 
is impossible. Deciding between these alternatives requires great numerical precision. 

The discontinuity of the effective action is sensitive to the definition of the effective 
degrees of freedom, just as singularities in the renormalization group may be created or 
eliminated by different choices of the block-spin transformation. A more sophisticated def- 
inition of (Tn, perhaps using a Kadanoff kernel to associate it with the smeared L^, may 
restore continuity, even without marked increase in the number of interaction terms. Note 
also that a reduction of the gauge theory to Z(3) spins in p resulted in an action that is con- 
tinuous across the phase transition. Perhaps a more complex effective spin, combining Ising 
with Z{3), will yield an effective action that offers both continuity and a local description 
of confinement physics. 
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FIGURES 



FIG. 1. Distribution of the Wilson line L^, averaged over m x m x m cubes, in the complex 
plane. Left: (3 = 5.0 (disordered phase). Right: /3 = 5.2 (ordered phase). Top to bottom: m = 1 
(single site), m = 2, m = 3. The lattice size is 2 x 8^. 

FIG. 2. Ln distributions for m = 2, as in Fig. 0, with circle |LnP = ~ ^-^ superimposed. 
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TABLES 




TART T? T Ot-iot'q+1-ivc tn' 


* appearing in the truncated effective Ising action. 


single spin 




n 


nearest neighbor 




n ij. 


next-nearest neighbor 




n n<i> 


3-spin bent 




n ii<u 


3-spin straight 




n n 


3'^'* neighbor 




n 


4th neighbor 




n /i 


TABLE II. Couphngs (3a for the (tentative) effective Ising action on a 16^ lattice, for the 
ordered and disordered phases at /? = 5.091. 


a 


ordered 


disordered 


1 
2 
3 
4 
5 
6 
7 


0.054(4) 
-0.455(1) 
-0.052(1) 
-0.0056(6) 
0.0063(9) 
0.044(1) 
0.152(1) 


-0.135(16) 
-0.390(8) 
-0.026(2) 
0.021(1) 
0.022(2) 
0.045(4) 
0.132(3) 
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TABLE III. Averages of C° in ordered and disordered phases of the gauge theory at /? = 5.091, 
compared with results of Ising Monte Carlo for the couplings listed in Table ||. Averages are 
normalized to 1. 



a 


ordered phase 




disordered phase 




gauge theory 


ismg ML/ 


gauge theory 


Ising MC 


1 


0.185(7) 


-0.849(2) 


-0.877(1) 


-0.696(6) 


2 


0.500(1) 


0.812(2) 


0.837(1) 


0.708(4) 


3 


0.353(1) 


0.775(2) 


0.804(1) 


0.626(5) 


4 


0.137(6) 


-0.738(3) 


-0.779(1) 


-0.595(6) 


5 


0.141(6) 


-0.742(2) 


-0.781(1) 


-0.602(6) 


6 


0.269(1) 


0.756(3) 


0.789(1) 


0.583(6) 


7 


0.206(2) 


0.743(3) 


0.781(1) 


0.552(7) 



TABLE IV. Couplings Pa for the effective Ising action after 


decimation to an 


8'^ lattice, for 


several gauge couplings 


surrounding 


the phase transition. 






a (5 = 5.08 


(5 = 5.091 


(3 = 5.091 (3 


= 5.1 


/3 = 5.2 




(disordered) 


(ordered) 






1 0.14(6) 


0.02(5) 


-0.027(2) 


-0.054(2) 


-0.26(3) 


2 -0.21(2) 


-0.24(2) 


-0.132(1) 


-0.131(1) 


-0.12(1) 


3 -0.043(3) 


-0.050(3) 


-0.023(1) 


-0.021(1) 


-0.025(3) 


4 -0.011(2) 


-0.015(2) 


0.0043(4) 


0.0028(4) 


-0.0014(16) 


5 0.003(5) 


-0.001(7) 


-0.003(2) 


-0.0003(9) 


-0.004(4) 


6 -0.002(3) 


-0.002(4) 


-0.006(1) 


-0.005(1) 


-0.006(2) 


7 -0.003(2) 


0.003(7) 


0.0016(7) 


0.005(1) 


0.004(3) 



TABLE V. Comparison of operator averages (0°) in ordered and disordered phases of the 
gauge theory at /3 = 5.091, compared with the effective Ising actions S'coid and S'hot simulated 
directly. 



a 


gauge theory 


"S'coid 


•S'hot 




gauge theory 




(ordered) 




(a)>0 


(a)<0 


(disordered) 


1 


0.188(7) 


0.154(2) 


0.9866(1) 


-0.8799(1) 


-0.878(1) 


2 


0.206(2) 


0.203(1) 


0.9739(1) 


0.7849(2) 


0.782(2) 


3 


0.123(2) 


0.119(1) 


0.9734(1) 


0.7772(2) 


0.774(2) 


4 


0.060(3) 


0.049(1) 


0.9615(1) 


-0.7006(2) 


-0.698(2) 


5 


0.063(3) 


0.053(1) 


0.9614(1) 


-0.6994(2) 


-0.697(2) 


6 


0.092(2) 


0.086(1) 


0.9733(1) 


0.7751(1) 


0.772(2) 


7 


0.080(2) 


0.073(1) 


0.9733(1) 


0.7745(2) 


0.772(2) 
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